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Abstract 

This paper shows the Fermi-Dirac Integrals Fi (rj) and Fs (rj) ex- 
pressed in terms of Riemann and Hurwitz Zeta functions. This is done by 
defining an auxihar function that permits rewrite the Fermi-Dirac integral 
in terms of simpler and known integrals resulting in the Zeta functions 
mentioned. The approach used here evades the use of iterative methods for 
the integrals and presents a clever procedure for around < 5, comple- 
menting Sommerfeld lemma, one that can be generalized for any integer 
k in Fk, (ri) in the refereed interval. 

1 Introduction. 

The motivation of this paper is to find expressions in order to calculate The 
Fermi-Dirac integrals specifically for 77 < 5, values that so far had been cal- 
culated using iterative methods as the trapezoidal method, [S]. This integrals 
appears in partial degenerate stars, |21, thermal conductions by electrons [1], 
and condensed matter physics. The Fermi-Dirac integral Fk{ri) is defined by 



nW=/ TTTT^d^ (1) 



Where 77 — ^r^, is the Boltzmann constant and ^ the chemical potential. 

In this paper we focus on the Fermi-Dirac integral Fi (rj) , and the procedure 
and result can be extended to the Fs (77) and generalized for any integer k in 
equation ([1]). This procedure involves an auxiliar function that relates Fi{r]) 
with ^"1(0), and based on their (graphical) relation will permit to write the 
auxiliar function in simple terms and make the integral Fi (77) easier to calculate 
obtaining Riemann and Hurwitz Zeta functions. That is, with almost no need 
for iterative methods. 



2 Fermi-Dirac Fi{r]) and the auxiliar function 

As we mentioned, we work on the Fermi-Dirac integral Fi (77) , defined by 



1 



Fi (??) = / ^ % (2) 



Actually, we focus on the term defined here by fi{(,, rf): 



m^v) = YT^ (3) 



and define the function /i(£): 



We define a function, denoted as /(£,??), that is the ratio of both functions: 



fi (rj) I + 

this function, graphically, behaves as a — be-"^ for ry < 5. That IS, we assume 
that we can model and define the function /(£,??) as: 



f{^,rj)=a-be-^ (6) 

where a, b, and c are constants. The a and b can be found considering some 
aspects of equation ([5]): 

/(0,77) = ^ 

and 

lim f{0,v) = e^ 
this implies that a = and h — • This means that 



= (7) 
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This leaves the issue of finding c. We take the following approach, perhaps 
the only one in this procedure that needs iterative methods: we find the max- 
imun of equation That is, after the algebra, we must solve for ^ 

in: 



+ e^{l-2^) =0 (8) 

and call this value f,„, that depends on 77. Using equations (O and (O, and 
substituting we can solve for c: 



Ln 

Cm 



e") - 1 



(9) 



It can be shown that calculating the two maximums, of equation ^ and 
/(^) V) * /i (?) will lead basically to the same result. Having found the constants 
in equation ([6]), we can write equation ([2]) as follows: 



1 + ei-^ 



dC 



- 1 
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—fdi - 



e-'t + lJo 1 + e« 



where the first integral in the equality we know from Arfken,]!], can be 
written as: 



1 - 

where is the well known Riemann Zeta function. For the second integral, 
after expanding some terms: 



f 

Jo 



i-iy 



= r(i) [2-k(i + i^)-C(i + ic + i) 

In the last equality we recognize C{p, q) as the Hurwitz Zeta function, where 
we used the identity as appears in Williams & Nan-yue, [Q. With these results, 
we can finally write Fi (77) as: 



(l-2-*)c(l.i) 



e-*? + 1 



, 1 c+1 



(10) 



C 1 + 
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3 Some results 



Using equation (jlOp . we elaborate in Table 1, for some rj, a comparison with 
values taken from [2 , which were based on McDougall & Stoner, Phil. Trans. 
Roy. Soc, 237:67(1938), including some values for positive 77 < 1, used in elec- 
tron conduction opacity, [3], showing how reliable and potentially applicable 
this method is. 



Table 1: Fi{ri) values comparation 



V 


Fi 

eq. cni) 


Fi 
2 

Clayton 


error (%) 


-4 


0.0161393 


0.0161277 


0.0715548 


-3 


0.0434453 


0.0433664 


0.181969 


-2 


0.11506 


0.114588 


0.411671 


-1 


0.292405 


0.290501 


0.655385 





0.678094 


0.678094 


6.95512E-11 


.1 


0.732034 


.733403 


0.18664 


.5 


0.977945 


.990209 


1.23858 


1 


1.35129 


1.39638 


3.22903 


2 


2.30003 


2.50246 


8.08906 


3 


3.58315 


3.97699 


9.90297 


4 


5.5495 


5.77073 


3.83358 


5 


8.99919 


7.83798 


14.8152 



4 Observations 

For T] = 1, the approach here presented is an improvement for the Trapezoidal 
scheme, [5], where their method gives an 6.4% error. The model according to 
equation ^ is reliable for 77 < 5, specially with negative 77. Yet, as the ratio 
of the functions in equation ^ shows a higher inflexion point, for r] > 5, the 
model needs an improvement. The reason for the limitation around this value 
appears when analyzing the equation (jlOp . where a simple examination shows 
a dependency of order , as it happens in Sommerfeld lemma. 

5 Conclusion 

The equation (jlOp gives an expression for calculating the Fermi-Dirac integrals 
Fi (77) in terms of the Riemann and Hurwitz Zeta functions. Using the same 
procedure we find expressions for -^3(7/) and Fk{ri), for general integer k, valid 
for 77 < 5. This approach complements that of Sommerfeld lemma, giving a 
nearly complete expression for Fermi-Dirac integrals. 
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^'1 



e'' - 1 
e-n + 1 



, 3 c + 
^+2 — 



(11) 



and for general k: 



i^.(77) = r 1 + 



)['•'( 



1-2-^ C 1 



e'' - 1 
e^'' + 1 



2-^1 



fc c + 
2'~ 



-C( 1+ 2>c+l 



(12) 



with the remark that ^m, equation ([5]), in order to find c (notice that c and 
/(^,?7) are k- independent), must be generahzed too: 



fce'' + e«(fc~20 = (13) 
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